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Abstract 

Caustic rings of dark matter were predicted to exist in the plane of the Galaxy at radii a„ ~ 
40 kpc/n for n = 1,2,3... The recently discovered Monoceros Ring of stars is located near the 
n = 2 caustic, prompting us to consider a possible connection between these two objects. We 
identify two processes through which the Monoceros Ring of stars may have formed. One process 
is the migration of gas to an angular velocity minimum at the caustic leading to enhanced star 
formation there. The other is the adiabatic deformation of star orbits as the caustic slowly grows 
in mass and radius. The second process predicts an order 100% enhancement of the density of disk 
stars at the location of the caustic ring. 

PACS numbers: 95.35.+d 
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I. INTRODUCTION 



The Monoceros Ring is an overdensity of stars in the plane of the Galaxy at a galacto- 
centric distance of approximately 20 kpc. It was discovered by the Sloan Digital Sky Survey 
fSDSS) ^ 3 and its existence was promptly confirmed by two independent collaborations 
[sl, 3]. The Monoceros Ring has been observed over 170 degrees in galactic longitude I in the 
Galactic anti-center direction (100° < I < 270°), and appears to be circular. Assuming it is a 
complete circle, the total mass in the Ring is estimated to be in the range 2 ■ 10^ — 5 ■ 10^ Mq 
in ref. j2[, and 2-10^ — 10^ Mq in ref. j^]- The scale height of the Ring stars in the direction 
perpendicular to the Galactic plane is estimated to be 1.6 ± 0.5 kpc in ref. jl], 0.75 ± 0.04 
kpc in ref. [sf and 1.3 ± 0.4 kpc in ref. 0]. The scale height in the direction parallel to the 
plane is also of order kpc. The stars in the Ring move with speed approximately 220 km/s 
in the direction of galactic rotation [El, 0] . Their velocity dispersion along the line of sight 
is small. It was estimated to be between 20 and 30 km/s in ref. 0, and 20 ± 4 km/s in ref. 

The most widely discussed interpretation of the Monoceros Ring is that it is a stream 



^9|. An alternative 
, [ll|. In this paper we 



of stars resulting from the tidal disruption of a Galactic satellite 
proposal is that the Ring is a manifestation of the Galactic warp [l( 
explore a different proposal altogether, namely that the Monoceros Ring of stars formed as 
a result of the gravitational forces exerted by the second caustic ring of dark matter in the 
Milky Way |12]. 

Caustic rings of dark matter had been predicted [l^] , prior to the discovery of the Mono- 
ceros Ring, to lie in the Galactic plane at radii given by the approximate law a„ ~ 40 kpc/n 
where n = 1, 2, 3 ... Since the Monoceros Ring is located near the second (n = 2) caustic 
ring of dark matter, it is natural to ask whether the former is a consequence of the latter. 
In our proposal, the position of the Monoceros Ring in the Galactic plane and its 20 kpc 
radius are immediately accounted for. 

Dark matter caustics are an unavoidable consequence of the fact that cold coUisionless 



dark matter (CDM) lies on a 3-dim. hypersurface in 6-dim. phase space [IJ, ll5|, [ly, [17 
The flow of dark matter in and out of the gravitational potential well of a galaxy necessarily 
produces caustics. These caustics are of two types, inner and outer. The inner caustics 
are rings when the angular momentum distribution of the dark matter is dominated by net 
overall rotation [l^ . The singularity structure of caustic rings was discussed in detail in ref. 
The ring radii were predicted in ref. 



15 



13l | using the self-similar infall model of galactic 



halo formation [l^, 20 1, generalized to allow the dark matter particles to have non-zero 
angular momentum 2l|] . Evidence for caustic rings of dark matter distributed according to 
the prediction of the self-similar infall model was found in the rotation curves of exterior 
galaxies [22] and the rotation curve of the Milky Way [l2| . 

Fig. [1] shows the transverse section of a caustic ring. Fig. 1(a) shows the flow of dark 



matter in the neighborhood of the caustic, whereas Fig. 1(b) shows the definition of its 
radius a and its transverse sizes p and q. As was mentioned already, the self-similar infall 
model predicts that the radius 02 of the second caustic ring of dark matter is approximately 
20 kpc in our galaxy. The transverse sizes p and q are not predicted by the self-similar infall 
model. However, for reasons explained in Section II, the expectation for p and q is that they 
are of order 1 kpc for n = 2. So the transverse sizes of the second caustic ring of dark matter 
are of order the transverse sizes of the Monoceros Ring. Moreover, for q = 1 kpc, the caustic 
ring mass enclosed within the triangular shape of Fig. 1(b) is approximately 6 ■ 10^ Mq (see 
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Section II). Thus the total mass of the n = 2 caustic ring of dark matter is of the same order 
of magnitude as the observed total mass of the Monoceros Ring. 

One might ask whether it is possible to interpret the Monoceros Ring as a caustic in a 
flow of stars. Indeed it is reasonable to assume that the flow of dark matter that forms the 
n = 2 caustic ring is accompanied by a flow of stars occupying the same 3-dim. hypersurface 
in phase space. Thus the proposal that the Monoceros Ring is a caustic in a flow of stars 
would explain equally well why the Monoceros Ring is in the Galactic plane and why its 
radius is 20 kpc. However, the proposal runs into difficulties. The first difficulty is that the 
self-similar infall model, which predicts the 20 kpc radius of the n = 2 caustic ring, also 
predicts that the matter in that caustic moves with a velocity of approximately 515 km/s 
in the direction of Galactic rotation (see Section II). This is inconsistent with the 220 km/s 
observed velocity of Monoceros Ring stars in the direction of galactic rotation. A second 
difficulty is that the matter in the caustic ring has an effective radial velocity dispersion of 
order 100 km/s, whereas the radial velocity dispersion of Monoceros Ring stars is only of 
order 20 km/s. Thus the proposal that the Monoceros Ring is a caustic of stars appears 
untenable. As mentioned already, the proposal we explore instead is that the Monoceros 
Ring is caused by the gravitational field of the n = 2 caustic ring of dark matter. 

We will find that there are two apparently viable mechanisms by which the n = 2 caustic 
ring of dark matter may cause the Monoceros Ring of stars. The first mechanism is enhanced 
star formation at the caustic radius because viscous forces on the gas in the neighborhood 
of the caustic drive it towards r = a. The second mechanism is the adiabatic deformation of 
the orbits of ordinary disk stars by the caustic ring of dark matter. We will show that the 
increase in star density at the location of the caustic as a result of this second mechanism 
is of order 100%. 

The paper is organized as follows: In Section II, we give a self-contained and hopefully 
pedagogical description of the properties of caustic rings of dark matter, including the pre- 
dictions of the self-similar infall model for the parameters that characterize the caustic rings, 
the observational evidence for caustic rings of dark matter, and the expected properties of 
the n = 2 ring. In Section III, we describe the gravitational field of a caustic ring of dark 
matter and the perturbation it causes to the galactic rotation curve. We then analyze two 
mechanisms by which the second caustic ring of dark matter may cause the Monoceros 
Ring: the migration of gas to a sharp minimum in the angular velocity at r = a implying 
enhanced star formation there, and the adiabatic deformation of ordinary disk star orbits 
by the slowly growing caustic. Section IV provides concluding remarks. Axial symmetry is 
assumed throughout unless stated otherwise. 



II. CAUSTIC RING PROPERTIES 



Cold coUisionless dark matter (CDM) particles lie on a 3-dim. hypersurface in phase 
space [13, HQ. We refer to this 3-dim. hypersurface as the "phase space sheet" . At the 
location of a galactic halo, the phase space sheet is wound up so as to cover physical space 
multiple times. This phase space structure implies that the velocity distribution of CDM 
particles is everywhere discrete and that there are surfaces in physical space, called caustics, 
where the density of dark matter particles is very large. Discrete fiows and caustics are a 



robust prediction of cold dark matter cosmology [17|. The reader may wish to consult ref. 



17( 1 for background information and a list of references. 
Galactic halos have outer caustics and inner caustics. 



The outer caustics are a set of 
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simple fold (A2) catastrophes located on topological spheres surrounding the galaxy at radii 
of order hundreds of kpc. Our focus, however, is on the inner caustics, which are much 
closer to the galactic center. The physical shape and catastrophe structure of inner caustics 



depend on the angular momentum distribution of CDM particles falling onto the halo [18[ . If 
that angular momentum distribution is dominated by net overall rotation, the inner caustics 
are a set of ring-like closed tubes, called "caustic rings", in or near the galactic plane. In 
cross-section each tube is a section of the elliptic umbilic (-D-4) catastrophe |15|. 

In this section we give a detailed desciption of caustic rings in the limit of axial symmetry 
and where their cross- sectional sizes, p and q, are much smaller than their radius a. Under 
these assumptions, the distribution of CDM particles in the vicinity of the caustic ring is 
determined in terms of a relatively small number of parameters, which we identify. Next, 
we summarize the predictions of the self-similar infall model for the caustic ring properties. 
We briefly review the evidence for caustic rings of dark matter, at the radii predicted by the 
self-similar infall model, in the Milky Way and in other isolated spiral galaxies. Finally, we 
list the expected properties of the second caustic ring of dark matter in our galaxy. 



A. Catastrophe structure 



The caustic ring singularity was analyzed in ref. [isl] which the reader may wish to consult 
for details. In the limit of axial symmetry and where the transverse sizes, p and q, of a caustic 
ring are much smaller than its radius a, the distribution of CDM in the vicinity of the caustic 
is given by the particle positions 



p(a,r) 



bar 



(1) 



We use cylindrical coordinates {z, p, 0) for position in physical space. Eqs. ([H) give particle 
positions at a particular time, say t = 0. The particles are labeled by parameters {a,T). 
a = ^ ~ 6 where 6 is the polar angle of the particle at the time of its last turnaround, r is 
the time when the particle crosses the z = plane, t — t can be thought of as the age of 
the particle. The particles labeled {a, r) form a circle of radius p{a, r) at a height z{a, r) 
above the z = plane. 

b, a, u , Tn and s are constants characterizing the caustic ring. Each has physical meaning. 
See ref. 15|] for more precise descriptions than we give here, a is the radius of the caustic 
ring. I To I is of order the time a constituent particle spends in the caustic, b is of order 
the speed of the particles in the caustic, u is of order their centrifugal acceleration, s 
characterizes the a-dependence of specific angular momentum near the equator (a = 0). 

Fig. [1] describes the caustic ring cross-section. Fig. 1(a) plots r), z{a, r)) for contin- 



uous r, and discrete values of a. The lines in Fig. 1(a) are the trajectories of the particles 
forming the flow, except that positions are plotted as a function of age, whereas for ordinary 
trajectories position is plotted as a function of time. Let us call the lines of Fig. 
trajectories" . 



Fig. 



la 



la 



•age 



shows that particle density diverges on a closed line which has the 
shape of a isosceles triangle, but with cusps instead of angles. We call that shape a"tricusp" . 

It is the envelope 



la' 



The location of the tricusp is shown in Fig. 1(b) for the flow of Fig. 
of the age trajectories. There are four flows everywhere inside the tricusp and two flows 
everywhere outside. The caustic, i.e. the surface where the density diverges, lies at the 
boundary between the region with four flows and the region with two flows. 
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The physical space density is given by 



1 ^i^"^ dM , . cosa 

{aj{p,z),Tj{p,z)) 



P j=i dndr \D(a,T) 

where aj{p,z) and Tj{p,z), with j = 1 ... N{p,z), are the solutions of p(«,t) = p and 
z{a, t) = z. N{p, z) is the number of flows at position (p, z); thus, = 4 inside the tricusp 
and N = 2 outside. D{a, r) is the Jacobian determinant of the map (a, r) (p, z): 

D{a, t) ^ det i^^] = -b[urir - To) + sa'] . (3) 

= cosadadT mass falling in per unit solid angle and unit time. The tricusp 

perimeter is the locus of points (p(a, r), z{a, r)) for which D{a, r) = 0. We call p and q the 



sizes of the tricusp in the p and z directions respectively; see Fig. 1(b), They are given by 



1 2 V27 b 

P=-^ur, , q=^-^^P . (4) 



One of the cusps of the tricusp points away from the galactic center. 

Consider the 3-dim. space spanned by the two physical coordinates p and z plus the 
constant tq. tq can be positive or negative. In (p, z, tq) space, the tricusp is perpendicular to 
the To axis with one cusp pointing in the positive p direction. As tq varies from negative to 
positive values, the size of the tricusp, which varies as Tq , shrinks to zero and then increases 
again. The structure in the neighborhood of tq = is the (full) elliptic umbilic catastrophe 
(D^/i), one of the elementary catastrophes in three dimensions. Thus the tricusp, i.e. the 
cross-section of a caustic ring, is a section of the elliptic umbilic. 

If the 2;-axis is rescaled relative to the p-axis so as to make the tricusp equilateral, the 
tricusp has a symmetry consisting of rotations by multiples of ^ about the point of 
coordinates (pc, Zc) = (a+p/4, 0). This point may thus be called the center of the tricusp. It 



is indicated by a star in Fig. 1(b) , As an example of how the above formalism is to be used, 
let us derive a formula for the density dc at the center of the tricusp. Using Eqs. ([T]), one finds 

that the four flows there have parameter values (a, r) = (0,|ro), (0,|ro), {.+\f^,0) and 

(~\/lf ' '^^^ corresponding values of the Jacobian determinant are D = —^bp, ^bp, —\bp 
and —\bp. Inserting these in Eq. ([2]), one finds 

abp dVtdr 

In obtaining this result, we approximated p by a and cos a by 1 in Eq. ([2]), and neglected 
the a-dependence of These approximations are appropriate since we assume that p 

and q are small relative to a. For the mass per unit length enclosed within the tricusp, we 
find by numerical integration 

A = 0.6 pqdc . (6) 

Note that some of the mass associated with the caustic ring lies outside the tricusp. In 
particular, when p = g = 0, A = but there is still a caustic ring and an associated 
overdensity. In that limit, the tricusp has collapsed to a point, and the particle density 
diverges as the inverse distance to that point [15i] . 
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Finally, we describe the velocities of the particles that constitute the caustic ring. Let us 
write the velocity of particle labeled (a, r) as v{a, r) = v^{a, r)0 + fp(a, r)/) + f^(a, t)z. For 
p,q a, all the particles in the caustic have the same speed v because they have come to the 
caustic from a small region (near the equator) of the turnaround sphere, along neighboring 
trajectories. The main component of velocity is in the direction: Vfj^ ~ v. The speed v is 
related to the centrifugal acceleration u by 

u^- . (7) 

In the case of a stationary flow, the velocity components in the p and z directions are 

dp dz 
Vp = -— = -u{t - To) , = -— = -ha . (8) 

Here we use the fact that, in case of stationary flow, the particle positions are func- 
tions only of their age t — t. For the four flows within the tricusp, Vp and Vz are 



of order utq = vy-^. For the four flows at the center, for example, {Vp,Vz) = 

(-|Mro,0), (+iMro,0), (utq, -I^utq), (utq, +1^uto), where ( = f = 

Caustic rings grow in mass and radius on cosmological time scales, implying that the flow 
is not exactly stationary. However, provided one has information on the manner in which 
the growth occurs, the time dependence is easily included. For example, if the growth in 
size is tantamount to an increase in the caustic radius at the rate d, the velocities are still 
as described in the previous paragraph except Vp is replaced with Vp + a. If the growth is 

tantamount to expansion in all directions by a scale factor R{t), Vp is replaced by Vp + ^p 
and Vz by Vz + ^z. 

Thus far, we have given a detailed description of a caustic ring in the limit of axial 
symmetry and where the transverse sizes p and q of the ring are much smaller than its 
radius a. The description is in terms of a small number of parameters: a, b, u, tq, s and 
-4M-. We now turn to the self-similar model of galactic halo formation to obtain estimates 
of these parameters for the caustic rings in actual halos. 



B. Predictions of the self-similar infall model 



The self-similar model of galactic halo formation 1^, 2^, 21 1 assumes that the entire halo 
phase space distribution is unchanged in time except for a rescaling of all lengths by a scale 



factor R{t), and all velocities by — where t is cosmic time. Physical space densities scale 
as ^. A spherically symmetric overdensity in an Einstein-de Sitter universe (^matter = 1) 
has self-similar evolution provided its initial profile is a power law 19, 2^ 



Mo 



(9) 



where Mj is the mass interior to initial radius rj, 5 Mi is the corresponding extra mass, and 
e is a parameter with the a-priori range < e < 1. The scale factor R{t) is proportional to 



2 

t3- 



The rotation curve is flat at small r provided < e < ^ [19[ . In an average sense, e is 



related to the slope of the power spectrum of density perturbations on galactic scales [23 
The standard CDM power spectrum implies e is in the range 0.2 to 0.35 21 
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The accelerated expansion [2J] of the universe is not consistent with Einstein-de Sitter 
cosmology, nor therefore with strict self-similarity. This is not a serious shortcoming of 
the self-similar model, however, because galactic halos were formed for the most part long 
before the onset of accelerated expansion, when the universe was accurately described by 
Einstein-de Sitter cosmology. 

The original spherically symmetric self-similar model [l^, 2^ assumes radial orbits for 
all the particles. Whereas this approximation is reasonable when describing the outer parts 
of a galactic halo, it is inadequate for the inner parts. Specifically, in the approximation of 
radial orbits, all particles pass through the galactic center each time they fall in and out of 
the galaxy, causing the density to diverge at the center as ^. The halo contribution to the 
rotation curve then approaches a constant as r — > 0. In actual galaxies, the central parts 
are dominated by baryons, and the halo contribution to the rotation curve goes to zero at 
the center. 



However, the self-similar model can be generalized [21| to allow the dark matter particles 
to have angular momentum. Self-similarity is maintained provided the specific angular 
momentum distribution on the turnaround sphere at time t is of the form: 

i{f,t)=j{f)^ , (10) 

where R{t) is the turnaround radius and j(f) is a dimensionless and time-independent an- 
gular momentum distribution. When angular momentum is included, particle orbits avoid 
the galactic center, the inner parts of the halo are depleted, and the halo contribution to the 
rotation curve goes to zero at r = 0, as it should. In the Milky Way, approximately half of 
the rotation velocity squared at our location is due to dark matter, and half is due to bary- 
onic matter. This determines the average magnitude j of the rescaled angular momentum 



distribution j(f) to be of order j ~ 0.2 for the Milky Way halo [21 . 

The caustic rings of dark matter occur where the particles with the most angular momen- 
tum are at their distance of closest approach to the galactic center 15, 3|. In the spherically 



symmetric model, the orbits are radial and hence the distance of closest approach vanishes 
for all particles. In that case all the caustic rings collapse to a single caustic point at the 
center of spherical symmetry. In the more realistic self-similar model with angular momen- 
tum, the caustic rings have finite radii. The radii are predicted in terms of the maximum 
value jinax of the dimensionless angular momentum distribution j(f) 



{a„: 71 = 1,2,3,4,5...}^ / ^^rot , ™ j 

^ , , , , J ^ ^220 km/s/ V0.25/ ^ ^ 

where v^ot is the rotation velocity of the galaxy. Eq. (ITTl) is for the particular case e = 0.3. 
However, the a„ oc ^ approximate behaviour holds for all e in the relevant range 0.2 < e < 
0.35, so that a change in e can be compensated for, as far as the a„ values are concerned, 
by a change in jmax- We will assume e = 0.3 henceforth. The self-similar model does not 
predict jmax and allows a different jmax value for each galactic halo. 

The relationship between j and jmax depends, of course, on the j(f) distribution. Let us 
assume the simplest distribution consistent with net overall rotation and axial symmetry: 
j(^) = jmax0cosa, i.e. that the turnaround sphere is initially rigidly rotating. In that case, 
j = f jmax- Thus, for the Milky Way halo, the estimate j ~ 0.2 implies jmax ~ 0.25. 

In summary so far, the self similar infall model with axial symmetry and net overall 
rotation predicts caustic rings in the galactic plane at the radii specified in Eq. (fTTj) . It also 
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predicts the speeds f„ of the dark matter particles forming the caustic rings and the infall 

dM 



rates ■^^\n- The f„ are approximately n-independent: 



The infall rates are given by 



515 ^ { , ) ■ (12) 

s V220 km/sy ^ ' 

dM 1 „ vl. 



dVtdr V 



with 

{/„:n= 1,2,3,4,5,...} ~ (11, 4.6, 2.9, 2.1, 1.7, ...) x lO'^ . (14) 

Eqs. f ll3l) and (HM provide the prefactor which appears in the formulas for the density 
[Eqs. ([2]) and ([5])], up to the ratio |, which is of order one but which may differ from one by 
a factor two or so. The transverse sizes, p and q, of the caustic rings and the ratio | depend 
on relatively more subtle properties of the velocity distribution at last turnaround , and 
are not predicted by the self-similar infall model. 

Observational evidence has been found in support of Eq. (fTT!) in the Milky Way and other 
spiral galaxies. We briefly describe this evidence now. 



C. Summary of observational evidence 

Caustic rings of dark matter in or near the galactic plane cause bumps in the galactic 
rotation curve. In ref. j22| a set of 32 extended and well-measured rotation curves was 
analyzed to test the hypothesis that some of their bumps are caused by caustic rings of dark 
matter at the radii given by Eq. fill I) . For each rotation curve, the radial coordinate r was 
rescaled according to 

/220 km/s\ , , 

r ^ f = r '— , 15) 

V ^'rot / 

where Wrot is the rotation velocity implied by the curve. The rotation curves were then co- 
added. The combined rotation curve shows peaks at f ~ 20 and 40 kpc, with a significance 
of 2.6 a and 3.0 a respectively. The result suggests not only the existence of caustic rings 
of dark matter at the radii given by Eq. ffTTl) . but also that the distribution of jmax values is 
peaked near jmax — 0.27. 

The Milky Way rotation curve has a series of sharp rises at radii which agree at the 3% 
level [12'] with the caustic ring radii given in Eq. flTT]) . The inner North Galactic rotation 
curve of ref. [25|] has ten rises between r = 3 and 8.5 kpc, which may be identified with 
caustic rings n = 5, 6 ... 14. The rises are "sharp" in the sense that they start and end 
with discontinuities (kinks) in the slope of the rotation curve. Kinks in the rotation curve 
are predicted by caustic rings because the dark matter density diverges at the caustic. The 
outer Milky Way rotation curve is much less well measured. Nonetheless it has a prominent 
rise near r = 13 kpc which may be identified with the n = 3 caustic ring. Finally, the 
IRAS map of the Galactic plane in the direction of galactic coordinates (/, 6) = (80°, 0°) 
has a triangular feature whose position and appearance is consistent with the imprint of the 
caustic ring of dark matter nearest to us (n = 5) upon the gas and dust in the disk. 
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D. Expectations for the n = 2 caustic ring 



Eq. ffTTl) predicts that the n = 2 caustic ring has radius 02 — 20 kpc. Eq. f|T2|) predicts that 
the speed of the particles constituting the ring is approximately 515 km/s. The transverse 
sizes P2 and q2 are not predicted. However, it is reasonable to expect that the n = 2 caustic 
ring has properties similar to caustic rings n = 3 and n = 5 to 14, for which we have the 
observational evidence mentioned in the previous subsection. The p values of those rings can 



be read off from the widths of the corresponding rises in the Milky Way rotation curve [12 
For those "observed" rings, one finds that p/a ranges from 0.015 to 0.1, with an average 
value of 0.05. We thus expect p2 to be of order 1 kpc, and very likely between 0.3 and 2 
kpc. The rises in the rotation curve do not inform us about the values of q. However, the 
triangular feature in the IRAS map provides values of both p and q for the n = 5 caustic 
ring. In that case | ~ 1.5. We therefore expect q2 to be of order 1 kpc as well. 

The prediction for the dark matter density at the central point of the tricusp is obtained 
by combining Eqs. ([5]), f|T3l) and ffUl) : 

_3 Mq fV2\ /kpc' 



Using Eq. ([6]), we obtain the predicted dark matter mass per unit length enclosed within 
the tricusp 

x^^n-io'^ r^) (^] . (17) 



kpc \62/ \kpc 

Since 62 is of order f 2, the total mass 27ra2A2 enclosed within the tricusp is of order 6-10* Mq. 
Finally, from Eq. ([8]) we learn that the transverse velocity components Vp and Vz of the 



particles within the tricusp are of order (515 km/s)y-^ ~ 160 km/s. 



III. EFFECT OF A CAUSTIC RING OF DARK MATTER UPON BARYONIC 
MATTER IN THE DISK 

In this section, we first describe the gravitational field of a caustic ring of dark matter 
and its effect on the Galactic rotation curve, assuming that the caustic ring lies in the plane 
of the disk. Next we identify two possible mechanisms by which a ring of stars may form 
due to the presence of a caustic ring of dark matter. The first mechanism is the migration 
of gas to the caustic ring radius, leading to enhanced star formation there. The second is 
the adiabatic deformation of star orbits as the caustic ring slowly grows in size. 



A. Gravitational field of a caustic ring 

As before, we assume that the transverse dimensions, p and q, of the caustic ring are 
small compared to its radius a. In that limit, when calculating the gravitational field g of 
the caustic ring at a distance of order p or q from the ring, we may neglect its curvature and 
pretend that it is a straight tube. Then 

Up, z) = -2G J dp'dz' dip', z') .Jf~'''"~/J.2 ■ (18) 
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Using Eq. ([2]), changing variables {p',z') neglecting the (a, r) dependence of 

^|^(«, r) over the size of the caustic, and approximating cos a ~ 1, we obtain 

^ 2G dM f {p- p{a,T),z - z{a,T)) 

Eqs. ([1]) provide the functions p(a, r) and 2;(a;, r). The integral on the RHS of Eq. f|T9|) may 
be evaluated numerically for all p and z. In the galactic plane 

with ^ = || (as before) and 

where X = For C = 1 oiie can do the integral by analytical methods, with the result 



/(1,X) = for X < 

= -i + for < X < 1 

= +^ for X > 1 . (22) 



A plot of /(1,X) is shown in Fig. 2(a) /(1,X) has an upward kink at X = and a 
downward kink at X = 1. The kinks are due to the divergent behaviour of the dark matter 
density in the z = plane at r = a and r = a + p. The shape of I{(,X) with ^ 7^ 1 is 
similar to that of /(1,X). The amount by which I{(,X) rises between X = and 

X = 1 depends on (. A plot of is given in ref. [1^]. is near one for ( of order 

one, and AJ(1) = 1. 

Substituting Eq. f|T3l) into Eq. fl20|) . we obtain: 

2 

?? 7; — n 

Ur,0) = -2f-^-IiC, )r . (23) 

r p 

2 

Since ^ is the gravitational field of the galaxy as a whole, and | and /(C, ^^) are both 
of order one whereas / is of order a few %, Eq. fl2^ implies that the caustic ring of dark 
matter causes only a small perturbation on the local gravitational field. 

Consider a smooth halo which produces a fiat rotation curve with rotation velocity frot- 
If we add to the halo a caustic ring of dark matter of radius a, the gravitational field near 
r = a would be ^(r) = —g{r)f with 

9{r) = '-f[l + 2fll{C,'-^)] . (24) 

However, this is not what we want to do. Indeed, the caustic ring of dark matter is not an 
addition to the smooth halo. Rather, it is the outcome of redistributing dark matter already 
present in the halo. One can see that it is incorrect merely to add caustic rings of dark 
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matter to a smooth halo by noting that this would cause the rotation curve to rise. The 
rotation velocity squared would rise by the relative amount 2fnf^AI{(n) at the nth ring, 
for each ring. Instead, the redistribution of dark matter into caustic rings should be such 
that the rotation curve remains flat on average. The gravitational field in the vicinity of a 
caustic ring is therefore of the form 



9ir) = ^ f 1 + 2/^[/(C, ^) + Hir)] ] (25) 



where H{r) is a smooth function which steps down by the amount A/((^) as r increases from 
r <^ a to r a; for example, H{r) = —^AI{() tanh(^^). p' is expected to be of order a. 
Let us define 

J{r-a)^l[IiC/-^) + Hir)] . (26) 
p 

The gravitational field in the galactic plane near a caustic ring of dark matter is then 

g{r) = ^[l + 2fJ{r-a)] . (27) 



Fig. 2(b) shows what the function J(r — a) looks like qualitatively. 



B. Effect on gas 

Gas in circular orbits in the galactic plane moves with angular velocity 



^ir) = f-^ . (28) 

Unless Q is r-independent (rigid rotation) there is shear in the velocity field. The viscous 
forces which result from this shear cause radial motion of the gas, as discussed by Lynden- 
Bell and Pringle [2^. For the sake of completeness, we repeat here relevant considerations 
from ref. 28l . 



Consider a cylinder of radius r. The viscous force per unit area across the cylinder is 
pvr^ where p is the gas density and u its viscosity. Hence, the viscous torque across the 
cylinder is 

T = lixr / az (ovr——) r = ZTCuar —— (291 
J dr dr 

where a = J dz p is the mass per unit surface of the gas. Consider an annulus of width 6r. 

Its mass is 6m = 27Tr5ra. The torque on the annulus is 

r(r + Sr) - r(r) = Sr-^{27riyar^^) = ^SL (30) 

dr dr dt 

where 5L = 5mVl{r)r'^ is its angular momentum, and ^ = §1 + "^r^ is the time derivative 
"following the motion". We have = 0. Also, ^[f2(r)r^] = because the gravitational 
field is determined by the mass distribution of the Galaxy as a whole and only negligibly 
altered by the motion of the gas. We thus obtain 



1 fMr^f) 
ar iinr^) 



(31) 
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for the radial velocity of the gas. The gas obeys the continuity equation 



| + ;^(™-)-o . (32) 

Eq. fl3Tl) neglects the back-reaction pressure from the accumulation of gas at particular radii. 

Let us assume that the gas is distributed uniformly = ^ = 0) to start with, and use 
Eq. fl3Tl) to determine in which direction it is driven by the viscous forces. In a smooth halo 
with flat rotation curve, fl = ^ and hence Vr = If a caustic ring of dark matter lies in 
the disk, we have instead in the neighborhood of the ring 

Q{r) = ^[1 + 2/J(r - a)]^ = ^[1 + /J(r - a) + 0(f)] . (33) 

Inserting Eq. (l33l) into Eq. (13T|) . setting ^ = ^ = 0, and neglecting terms of order p, we 
obtain 



u l + fj{r-a)-rff-r 'f^ 
r l + /J(r-a)+r/^,^ 



To obtain the qualitative behaviour, we neglect /J(r — a), t/^ and r^/^^ vs. 1, because 



such terms are all of order f. Then 

V 1- fr^ - fr'^f^ 

dl 



r i + fr 



(35) 



dr 



For r < a and r>a + p, ^ = and hence = — ^ as in the absence of caustic. For 

a < r < a + p, 

1 - + 

(36) 



where X = as before. Since = 0(1) whereas ^ = 0(10), the numerator on the RHS 
of Eq. (1361) is always dominated by the last term. Hence 



V fa' 



zp 



(37) 



Eq. (!37|) shows that, within the tricusp, the gas velocity is also inward but much larger, by 
a factor three or more, than outside the tricusp. It increases rapidly as r ^ a+. Therefore 
the viscous forces tend to drive the gas towards r = a. This happens on gas dynamic time 
scales which are much shorter than the cosmological time scales over which the caustic rings 
migrate. The accumulation of gas at the caustic ring radius may lead to enhanced star 
formation there. Such processes may be at work in the case of the Monoceros Ring. To 
the extent that such processes dominate the formation of the Monoceros Ring, the stars in 



the Ring should be younger than average. Ref. [26| presents evidence that Monoceros Ring 



stars are on average bluer and therefore younger than ordinary disk stars. 



C. Effect on star orbits 

In this subsection, we consider the adiabatic deformation of disk star orbits by the slowly 
growing caustic ring of dark matter. 
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1. Collision rate 



Adiabatic deformation of star orbits presupposes that the star colhsion rate is smalL The 
time scale over which star orbits are significantly modified through gravitational scattering 



with other stars in a population is [27 



trclax ^ 0.3 , . . 38 

where n is the density of stars in the population, m their typical mass, a their velocity 
dispersion, and InA ~ 20. For the stars in the Monoceros Ring, treiax ~ 10^^ yr if we set 
m ~ Mq, (T ~ 20 km/s, and n ~ lO^^/pc'^. Since treiax is much greater that the age of the 
universe, we are justified in neglecting collisions among stars in the Ring. 



2. Orbit stability 

In a gravitational field ^(r) = —g{r)f, the angular frequency squared of small radial 
oscillations about a circular orbit of radius r is 

-'W = ^|:P^^W) • (39) 

The orbit is stable ii > 0. In the neighborhood of a caustic ring of dark matter, where 
the gravitational field is as given in Eq. (127|) . we have 



a;^(r)=2(^)'(^l + 2/J(r-a)+r/^) . (40) 



The sum of the first two terms in the parentheses on the RHS of Eq. (1401) is of order one. 
The third term is 

V (dl , dH\ , 

^ is everywhere positive, whereas ^ has the qualitative form 

dr ^ 2p'cosh^(^) ' ^ ' 

Since we expect p' to be of order a and at any rate much larger than fa ~ ^a, we find 
that circular orbits are stable everywhere in the neighborhood of a circular caustic ring of 
dark matter. The conclusion is valid only if the caustic ring is circular, because lack of axial 
symmetry may drive an instability through the phenomenon of Lindblad resonance, as we 
now discuss. 

A particle on an orbit of radius r is subjected to a radial time-dependent, but periodic, 
gravitational force from a non-circular caustic. The period T = ^ ~ equals the 

time to go around the Galaxy once. The equation of motion for the radial coordinate of 
the particle is that of a harmonic oscillator of proper frequency uj{r) driven by a periodic 
external force with frequencies mfi(r) ~ m— , where m = 1, 2, 3 ... . Resonance occurs 
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when uj{r) = mQ{r). For the corresponding radii, circular orbits are unstable. Using 
Eqs. ( !33|) and (140|) . we find the instability condition 

rf- = (— - 1) [1 + 2/J(r - a)] (43) 

in the neighborhood of a caustic ring. Since 2/J(r) ^ 1, and 

dJ dl 1 1 



the instability occurs at radii 

rm^a + f— ^ (45) 
p (m^ — 2)^ 

for m = 2, 3, 4... There is no resonance for m = 1. For a = 20 kpc and / = 0.046, the 
instabilities occur at 

{r^ - a : m = 2, 3, ...} = (210, 17, ...) pc ^ (46) 

Let us caution that these estimates take account only of the gravitational force exerted 
by the caustic ring of dark matter itself, and neglect the gravity of the baryons that have 
aggregated at the caustic. 



3. Circular orbits 

In the absence of the caustic, the effective potential for radial motion is 

Vesfiir) = vlJnr + — (47) 

where I is the specific angular momentum of the star. When the caustic is present, the 
effective potential is Vcs{r) = Vcs,o{r) + Vc{r) where 

V,ir) = 2fvl,lyJir-a) . (48) 

Vc{r) is plotted in Fig. 2c for 

J(r -a)= /(I, ^) - \ tanh (^-^] (49) 
p 2 \ P J 

with a = 20 kpc, p = 1 kpc, and p = 5 kpc. Fig. 2c illustrates the fact that the effective 
potential is smooth even though its second derivative diverges at r = a and r = a + p. 

The caustic ring radius increases with time. According to the self-similar infall model, 
a oc t3 + 97. Consider a particle of specific angular momentum /. In the absence of the caustic, 
it is on a circular orbit of radius r given by 

f = r^g{r) = r^vl, , (50) 
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or it is oscillating about a circular orbit of that radius. In the presence of a caustic of radius 
a, the particle is on or oscillating about a circular orbit of radius r (r, a) given by 



/2 = r'3^'(r')=r'2<, [l + 2/ J(r' - a) 
Angular momentum conservation implies 

= r (r, a)^ 1 + 2/ J(r (r, a) — a) 



(51) 



(52) 



Fig. [3] shows r' (r, a) as a function of a. Each line in that figure corresponds to a different 
value of r. Let p(r) be the density of stars in the absence of the caustic, and p (r, a) their 
density in the presence of a caustic with radius a. Assuming that all stars are and remain 
on circular orbits, conservation of the number of stars implies 



r dr p {r ,a) = rdr p{r) 
with r'(r, a) given by Eq. (!52|) . One readily finds 



p'(r',a) = p(r) 



1 + 2/ J{r' - a) + fr' ^{r 

dr 



(53) 



(54) 



Assuming that the initial star density has no significant structure of its own, we have 



p'(r,a) = p 



l + 2fJ{r-a) + fr^{r) 



(55) 



where p is the initial density. The second term in brackets is an order 10 % modulation of 
the initial density. The third term is of order /, like the second term, but diverges when 
r — >■ a+ as I/a/^ — a. It has the same singularity structure at r ^ a+ as the dark matter 

density because g' (r) = ^ [1 + 2/ J(r)] satisfies Gauss' law. Eq. (155!) states that, in the 
limit where all stars are on circular orbits, the star density adopts at r = a+ the same 
divergent profile as the dark matter in the caustic surface there. The pile-up of circular 
orbits at r = a_|_ is clearly seen in Fig. [31 



4- Non circular orbits 

The observed radial velocity dispersion Av of stars in the Monoceros Ring is of order 20 
km/s. This implies that the stars do not move on circular orbits, but oscillate in the radial 
direction with typical amplitude 

km 1 km 1 20 kpc , , ^ 

Ar ~ 20 - 20 ^ — = 1.3 kpc . (56) 

S UJ S V2 frot 

The density profile of stars in the neighborhood of a caustic ring is therefore the profile of 
Eq. fl5^ averaged over the length scale Ar. The sharp features at r = a and r = a + p are 
then washed out. However, there will be an average relative overdensity within the tricusp 



p -p 




af , kpc 
~ — ~ 1 (57) 

a<r<a+p P P 
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provided p is not much less than Ar. This is the 100% average overdensity of ordinary disk 
stars at the n = 2 caustic ring of dark matter, announced earher. 

For the proposed interpretation of the Monoceros Ring, the radial velocities acquired by 
stars as a result of the passing of a caustic ring should be less than the measured radial 
velocity dispersion of stars in the Monoceros Ring. We now verify that this is the case. 

The radial motion of a star in a near circular orbit in the neighborhood of a caustic ring 
is that of a harmonic oscillator whose equilibrium position is r (r(/), a{t)) and whose proper 
frequency u is given by Eq. fHOl) evaluated at r . a(t) is a smooth function of time. However, 

as shown in Fig. [3l r' as a function of a has kinks at a = r' and a = r' — p. The velocity ^ 
of the equilibrium position changes abruptly when it passes the caustic at those locations. 
A star in circular orbit in the z = plane first encounters the caustic at r = a + p. For 
illustrative purposes, we consider the case a = 20 kpc, to = 13.7 Gyr, e = 0.3, p = 1 kpc, 
V = b and / = 0.046. The caustic moves with outward speed 



da 
'dt 



2 

9^ 



a 

To 



2.0 km/s 



(58) 



When r' = a + p, the velocity of the equilibrium position changes abruptly by the amount 
{S is infinitesimally small) 



A 



dr 
H 



dr 
da 
af 



(r = a + p 
da 



6^ - ^ (r' = a + p + 5 



da 



0.6 km/s 



(59) 



2p + af dt 

Here and below we are neglecting terms of order /, p/a, p/p' and af /p' versus terms of order 
one. The proper oscillation frequency after the caustic has passed is 

1 



uj{r = a + p — 6) 




5 X 10'^ years 

Thus a star which is in a circular orbit in the 2; = plane before the r' 
passes by, will be oscillating in the radial direction with initial amplitude 



A 



1 A dr 
-A — 

uj dt 



30 pc 



(60) 

a + p caustic 

(61) 



after that caustic passes by. Stars not in the 2; = plane do not go through the caustic cusp 
at r = a + and presumably make smoother transitions through the caustic than the stars 
in the 2 = plane. So, Eq. fl6Tl) provides an upper limit for the 2; 7^ stars. Stars receive 
a second jolt when the r = a caustic surface passes by. The analogous quantities for that 
transition are : 



A dr 

^ -r 

dt 

uj{r = a_) 



dr , 
(r : 

da 

a 



dr 
da 
1 



da 
It 



da 



—2 km/s 



A 



A 



dr 
H 



1 



6 X 10'^ year 
130 pc 



(62) 



After the equilibrium position r has passed by the caustic surface at r = a, r decreases 
very quickly, implying a depletion on the r < a side of the caustic. This depletion is evident 
in Fig. [31 
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IV. CONCLUSIONS 



The flow of cold collisionless dark matter in and out of the gravitational potential of 
a galaxy necessarily forms inner and outer caustics. The inner caustics are rings in the 
galactic plane provided the angular momentum distribution of the infalling dark matter is 
characterized by net overall rotation. Assuming self-similar infall, the radii of the caustic 
rings of dark matter in our galaxy were predicted to be 40 kpc/n with n = 1, 2, 3 ... Because 
the Monoceros Ring of stars is located near the second caustic ring of dark matter we looked 
for processes by which the latter may cause the former. 

We have identified two such processes. The first is the flow of gas in the disk towards the 
sharp angular velocity minimum located at the caustic ring radius, increasing the rate of 
star formation there. To the extent that this process is responsible for the formation of the 
Monoceros Ring, the Ring stars are predicted to be bluer than average. The second process 
is the adiabatic deformation of star orbits in the neighborhood of the caustic ring. As the 
spatial dependence of the gravitational field of a caustic ring is known, it is straightforward 
to obtain the map of initial to final orbits for disk stars. The resulting enhancement of 
disk star density at the location of the second caustic ring is of order 100%. Because of 
uncertainties in the caustic parameters | and p, and in the velocity distribution of the disk 
stars, the strength of the enhancement can only be estimated within a factor of two or so. 

The self similar infall model of galactic halo formation is expected to describe the halos 
of all isolated spiral galaxies. The caustic rings of dark matter in exterior galaxies may also 
be revealed by the baryonic matter they attract. Our analyis is relevant to those cases as 
well. 

The existence of caustics has implications for most approaches to the detection of dark 
matter, including direct searches of WIMPs j29| and axions the gamma ray signal from 



WIMP annihilation [3l|], and gravitational lensing [32 
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(a) 




(b) 

FIG. 1: (a) Dark matter trajectories forming a caustic ring of dark matter, in p — z cross-section, 
(b) The envelope of the trajectories shown in (a). We refer to the shape shown in (b) as the 
"tricusp". The figure also indicates what is meant by the radius a of a caustic ring, and by its 
transverse sizes p and q. The star indicates the center of the tricusp, as defined in the text. 
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-o.i - 



-10 -5 ^ JC ^ 

(c) 

FIG. 2: (a) The function /(I, X) defined in Eq. 1^2), with X = (b) The function J{X) defined 
in Eq. for v = b, ( = I and H{X) = -itanh(f ). (c) The function Vc{X) defined in Eq. (j48]) . 
in units of for '^(^) as in (b). 
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16 18 20 22 24 



FIG. 3: The radii r'(r, a) of circular orbits in the presence of a caustic, as a function of the caustic 
radius a. r is the radius of the orbit in the absence of caustic. Each hne corresponds to a different 
value of r. Each line has an upward kink at r' = a + p and a downward kink at r' = a. As a 
increases with p fixed, r' first decreases from its initial value r until r' = a + p. For a < r' < a + p, 
i.e. when the orbit is within the tricusp, r' increases with a but not as fast as a, until a = r'. As 
a increases yet further r'{r,a) returns to the value r from which it started. Note that the radii of 
circular orbits pile up at r = a. Stars on approximately circular orbits oscillate about r'{r, a). 
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